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Abstract. We prove the convergence in a strong norm of a finite difference 
semi-discrete scheme approximating a coupled Schrodinger— KdV system on 
a bounded domain. This system models the interaction of short and long 
waves. Since the energy estimates available in the continuous case do not 
carry over to the discrete setting, we rely on a suitably truncated problem 
which we prove reduces to the original one. We present some numerical exam- 
ples to illustrate our convergence result. Nonlinear Schrodinger equation and 
Korteweg— de Vries equation and short wave long wave interaction and finite 
difference scheme 

AMS Subject Classification. 35Q55 and 35Q53 

1. Introduction 

In [5], D.J. Benney presents a general theory modeling the nonlinear interaction 
between short waves and long waves, deriving nonlinear differential systems de- 
scribing these interactions in various physical settings. The (complex- valued) short 
waves u{x,t), a; G K, f > 0, are described by a nonlinear Schrodinger equation and 
the (real- valued) long waves v{x,t) satisfy a quasilinear equation, eventually with 
a dispersive term. In the most general context, the interaction is described by the 
nonlinear system 

idtu + icidxU -\- dxxU — ctvu + j3\u\^u 
dtv + C2dxV + fid^v + vdxv^ = ^dx{\u\'^), 

where ci, C2, a, /3, 7, and v are real constants. 

In this paper, we are concerned with the numerical approximation of the solu- 
tions to the Cauchy problem for the system comprising the nonlinear Schrodinger 
equation coupled with a Korteweg-de Vries equation with Dirichlet boundary con- 
ditions on a bounded domain (0,i), 

(1.1a) idtU + dxxU — avu + (3\u\^u 

(1.1b) dtv + dlv + dx{v'') ^ -idx\u\^ 

(l.lc) M(a;, 0) = Mo(x), v{x,Q) ^ vq{x), xe(0, L), L > 0, 
(l.ld) M(0,i) = u(L,t) = 0, w(0,i) = w(L,t) = 0, ie[0,T], T > 0. 
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This kind of system arises in several fields of physics such as the study of resonant 
interactions, short and long capillary-gravity waves on water [8], an electron plasma 
interaction [II and a diatomic lattice system [13 . 

The Cauchy problem for the system ( 1.1 ) was initially studied on the whole line 
by M. Tsutsumi [T^], who proved that for initial data {uq,vo) in i/'"+^/^(M) x 
i?™(E), m — 1,2...., the problem is well-posed in the same space. After that, 
Bekiranov et al. [T] established local well-posedness for initial data in H^'iM) x 
s > 0, and more recently, Corcho and Linares ^ proved global well- 
posedness in the energy space x 

It is worth pointing out that in the above well-posedness results, uniqueness is 
obtained only in some subspace of C{[Q,T]; {M.) x H'-^/'^{m.)) (in P), (resp. a 
subspace of C([0, T]; ff^(M) x H'^{R)) in 4 ). Additionally, the techniques used in 
the papers^ i4j (which were introduced by Bourgain j3j and developed by Kenig, 
Ponce and Vega [21 UHl), do not seem to be appHcable to the generalized KdV 
equation (gKdV). In a very recent paper, Dias et al. [F, using energy methods, 
obtain a global solution in (iJ^(M))^ for a coupled Schrodinger-gKdV system. 

In this paper, we prove a convergence result for a semi-discrete finite difference 
approximation of the system (1.1) in the space (i?^(0,i))^. The energy methods 
used by M. Tsutsumi .V2^ to prove global existence of a solution fail in the discrete 
setting, so we propose a new approach: by an appropriate truncation of the qua- 
dratic function appearing in the equation (1.1b), we consider a related problem 
involving a gKdV equation. For each fixed value of the truncation parameter, we 
are able to prove the convergence of a numerical scheme toward the solutions of this 
auxiliary problem. These solutions, in turn, satisfy an energy inequality. Lastly, 
using this energy inequality, we are able to derive an L°° estimate independent of 
the truncation parameter, which implies that the truncated problem in fact reduces 
to the original one. 

In contrast to previous work, the proof of these stability estimates require that 
we work on a bounded subset (0, L) C M. From the viewpoint of the applications 
and numerical approximation, this presents no great loss in generality. 

Note also that our convergence proof does not rely on any previous existence 
results, and so constitutes a new existence proof for the Cauchy problem (1.1 1. 
Additionally, the present work represents, to the authors' best knowledge, the first 



numerical treatment of the system ( 1.1 1-( 1.1c 



An outline of the paper follows. After some notations and preliminaries, we 
state in Section [2] the main convergence result. Theorem |2.1[ In Section [3j we prove 
Theorem |2.1| and present Proposition 3.2 our main auxiliary result, dealing with 
the convergence of approximate solutions to a suitably truncated system, and an 
energy estimate. Its proof is the object of Section |4j Finally, in the last section 
of the paper we illustrate our convergence result with some numerical simulations 
and check its accuracy by testing it against some known exact solutions. 



1.1. Notations and preliminaries. Let us introduce the Banach spaces 
Xj^c = {z^ = {zq, zi, . . . , zj+i) e C'^+^ ; zo = zi = zj = zj+i = 0} 



with J e No and h — L/{J + 1). In a similar way, we define the real space Xjr. 
When no ambiguity arises, we represent either of these spaces by Xj. The scalar 
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product is given by 

j-i 

{z^,w^) = hzjWj, z^,w^ e Xj, 
i=2 

and the p- norms by 

i/p 

\\A\p,h= (y2h\zj\p) \ l<p<oo; |lz''Hoo= max \z,\, z'^eXj. 

\ ^ — ^ / 7— 2,...,J— 1 

To simplify the notation we write ||2;||p for the norm of z in both the continuous 
and the discrete case. We denote by H"'{0,L), H^{Q,L) and F-™(0,L) (m £ N) 
the usual Sobolev spaces. All the norms appearing in this paper are in (0,i), so 
for instance ||w||2 means |iu||L2(o_i). 

We will use the following notations for the finite difference operators. For z — 

D+z, _ - , D^z, - - , 



Doz, = ^-^±^^^^-{D^ + D_)z„ 



A^zj = D^D+Zj = D+D^zj = 



D^z, = DoD^D+z, = —■ 

For u = (uj), let us now introduce the piecewise linear interpolator, 

(1.2) P^u{x) = Uj + {x- Xj) ^''^^ _ , xe{xj,Xj+i), 

Xjj^l Xj 

and the piecewise constant interpolator, 

PoM(a:;) = Uj, x € {xj,Xj+i). 
Let u'^ G Xj. From (IT^ we derive 

(1.3) WP'lu'' -P'iu% < Ch\\D+u% 
for some C independent of h. As a consequence, we obtain 

(1.4) IIPXII2 < ||PX - Po'w'lb + WP'^u'^h < C\\u%. 

The following lemma establishes some inequalities which will be of use through- 
out. 

Lemma 1.1. Let (p — € Xj. Then 

(1-5) II0II.O < ^/2||0||^/^P±0||^/^ 

(1-6) < l\\D±ni < ^\\D±cf>h 

(1-7) ||<^||2 <q|0IU, 

where C is a constant independent of h. 
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Proof. The inequality ( 1.5 1 is a consequence of the Gaghardo-Nirenberg inequahty 
and d,,P'l 



iip?0iioo < V2\\pmy'\\d,p^,4,\\i/' < V2\ 



1/2| 



1/2 
2 I 



< 



The last 



while (1.6) is a consequence of the (continuous) inequality ||( 
inequality is an elementary consequence of the definition of the discrete norms. □ 

2. Statement of the main result 

We propose the following scmidiscrcte finite difference approximation to the 
Cauchy problem (1.1 1: 

(2.1a) idtu^ + A^u'' = (3\u^\^u'' + av^u^ 

(2.1b) dtv'' + D\'' + Doiv'^f = 7^ok''f 

(2.1c) u'^{0) = u'i, v\0)^v^, 

(2.1d) u''it),v''it)eXj, te[0,T]. 

The global existence proof of Tsutsumi [12] relies on energy methods which we 
cannot carry over to the finite difference framework. It turns out that in the semidis- 
crete case, the crux of our convergence argument relies on an a priori L°° bound. 



However, this bound is only available for a modified problem (see Proposition 3.2 



below). To deal with this difficulty, we use the fact that, under the right conditions, 
this problem reduces to the original one. 

Our main result establishes the convergence of the approximations (2.1 ) toward 
a global weak solution of the problem (1.1) in the space (_ff^(0,i))^. 

Theorem 2.1. Let a, /3,7 he such that 07 > 0, (uo,wo) G (-ffg(0, L))^ and let 
{u^^v^') e (C([0, T]; Xj))'^ he the global solutions of the discretized prohlem (2.1) 
with initial data {uq,Vq), such that PJ'mq 
/i — 5- 0. Then, up to a suhsequence, 

P^u'* ^ u, P^i;'' ^ V 1 
P'l'u'^ ^ u, V\v^ V 



uo and P^Wq — > vq in H^{0,L) as 



L°-{[0,T];H\0,L)), 
L°°{[0,T];L\0,L)), 



with 



(u,«)e (L-([0,T];i/i(0,L)))' 

n (C([0,r];i7i(0,L)) X C([0,T];L2(o,L))), 



T > 0, 



a weak solution of the Schrddinger-KdV system (1.1). 

3. Proof of Theorem 12.11 and stability estimates 



In this section we prove Theorem |2 . 1 1 along with the necessary stability estimates. 
We begin with the definition of an appropriate truncated problem. For each M > 1 
we define C°° functions f^^ , g^' satisfying 

if|w|<M, 
\vl if|i;|>M2 + l, 



and 



V, if \v\ < M, 
± C, if \v\ > 2M, 
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with < < v"^. Here, the constant C = C{M) is chosen to ensure the 

following property, 

(3.1) |(/*')'|oo + |.g*'|oo < C{M), |(/^)'U < 1. 

The functions /^^ , g*^ arc simply appropriate smooth truncations of the functions 



and V appearing in (1.1b I. We define also F^'{v) — f^^ {s) ds 



Now, we consider the auxiliary Cauchy problem, 
(3.2a) idtu + d^xU = (i\u\^u + ag'^' {v)u 
(3.2b) dtv + dlv + dj^'{v) = 7a.((/0'Mkn 
(3.2c) u{x,Q) ^ uq{x), v{x,0) = vq{x), xe(0, L), L > 0, 
(3.2d) u(0,i) = u(L,t) = 0, w(0,i) = v(L,t) = 0, ie[0,T], T > 0. 



and we propose the following semidiscrete finite difference approximation of (3.2): 
(3.3a) idtu'' + A'^u^' = p\u''\^u^ + ag^\v'')u^ 

(3.3b) dtv^ + + D^f^'^iv'^) = 7i?o((/')'(^^")|^^'f ) 

(3.3c) 
(3.3d) 



We will need some conservation laws of the auxiliary system (3.2). A first result 
is as follows: 

Proposition 3.1. Let I cM. be an interval and let {u,v) G {L^^{I; Hq{0, L)))'^ be 
a solution of the auxiliary system ( |3.2[ ). Then, for each t,s £ I, we have 

(3.4) M^\t) := \\u{t)h = \\u{s)h = M*^(,s) 



(3.5) Q^^{t):=a v\x,t) dx + 2-flin u{x,t)dxu{x,t) dx ^ Q'^\s). 
Jo Jo 



Proof. Although the estimates (3. 4), (3. 5 1 are formally easy to obtain, the rigorous 



justification of (3.5) is nontrivial and requires techniques from semigroup theory. 



First, since u{t) € Hq{0,L), from the equation (3.2a) we deduce that 

lm{idtu,u)H-ixHi +lm{9xxU,u)H-ixHi = 0. 

and so, 

d 



dt 



\uY{t) dx = 0, 



which gives (3.4 1 



For the proof of (3.5 ), we follow the ideas of Kato [71 Lemma 3.1]. We point out 
that —d^ is a skew-adjoint operator on L^(0, L) with domain H^(0,L) D Hq{0, L) 
and it generates a group of isometrics Uxit) on L^(0, L). Also, iA is a skew-adjoint 
operator on L^(0, L) with domain H^{0, L) n H^{0, L) and so it generates a group 
of isometrics Us{t) on L^{0,L). Now, we write the equations (3.2 1 in integral form. 



(3.6) 
(3.7) 



i{t) = Us{t - s)u{s) + / Us{t~r)m{r)dr, 



v{t) ^Unit- s)v{s) + / UK{t~r)n{r)dr, 
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with 



m(r) = —iag^ (v(r))u{r) — if3\u(r)\'^u{r), 
n{r)^-d,f^'{v{r))+^d4{g"'y{v{r))\u(r)\^). 



The formulas (|3^,([3J]) are verified in L'^{0, L) and, since d^m G L^^{I] ^^(O, L)), 
we have also 

(3.8) d,u{t) ^Us{t~ s)dMs) + J Us{t~ r)d^m{r) dr 

in i2(0,L). 



(3.8) 



Using the isometric property of Uxii) and Us{t), we easily deduce from (3.6|- 
{v{t),v{t)) = {v{s),v{s)) + 2 / {v{r),n{r)) dr, 



Im Im (-u(s), 9a;M(s)) + 2 J Im (m(r), 9^u(r)) dr. 

Since u(r),v{r) e Hq{0,L), we obtain from the expressions of m(r) and n{r) 

Im (^m{r) , dxu{r)) —— Ke a / g^\v{r))u{r)dxu{r) dx 



fSRe I \u{r)\ u{r)dxu{r)dx 

L 



It follows that 



{g^')'{v{r))dMr)Hr)\^dx. 
27lm (TO(r), — —a(v{r),n{r)'j , 



which implies the conclusion (3.5). This completes the proof of Proposition 3.1 □ 



The following result establishes the convergence of the approximations (3.3) to 
a global solution of the truncated problem (3.2), and the crucial energy estimate 
pl2| ). 

Proposition 3.2. Let a, (3, j be such that aj > 0. For each M > let {u''-'^'^ ,v'^-^'^) e 
{C{[0,T]; Xj))'^ be the solution of (3.3) with initial data {uq,Vq) such thafP^u^ — > 
Wo, PjwJ — > Wo in H^{0,L) as ft. — > 0. Then, up to a subsequence, 



(3.9) 

(3.10) 

with 

(3.11) 



h„.h,M 



P>' 



,M 



h„.h,M 



P> 



in L°°{[0,T];H\0,L)), 
m i°°([0,T];L2(0,L)), 



n (C([0,r];i7i(0,L)) X C([0,T];L2(o,L))), T > 0, 



a global weak solution of the truncated system (3.2). Moreover, the following energy 
estimate is valid, 

Z^\t) 



(3.12) 



- cxF^\v''\t)) + ^|u*'(t)|^} dx < £*^(0), 
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for all t e [o,r]. 

We po stpo ne the proof of Proposition |3.2| until Section |4] and proceed to prove 



Theorem 



2.1 



The goal is to prove an L°° bound for u 



M „,M 



independent of the 



truncation parameter AI, using the energy inequality ( 3.12[ ). Once this is achieved, 
it is clear from the definition of and g^^ that taking AI large enough yields a 
solution of the original problem, (1.1). 



Let us define Mq = M(0), Qo = Q(0) (see Q , ([3]5]) ) , and set 

\^\\d^v,\\l + \a^\\\v,\\,\\uo\\l+^-^ 



(3.13) £o - hmuoWl + ^||9.^;o||2 + l«7lll«o||2ho||4 + ^ll«o||^ + ^ll^ollt 



Observe that |£*^(0)| < £o for all AI > (sec ( [3l2| )). 

For simplicity, in what follows we omit the superscript AI from the solutions 



(u ,v ) of the system (3.2) obtained in Proposition 3.2 



First, note that the energy inequality (3.12) gives 



(3.14) 



< {E.o + \aj\Ht)h\\u{t)\ 



2 



Mt)\\t) 



Next, from (3.5) we have 
(3.15) 



\\vml<--{\Qo\ + 2Muoh\\dMt)h). 



a 



Let now m = min{|7|, |q!|/2}. Using again Gagliardo-Nirenberg and Young in- 
equalities, we deduce from ( |3.14 1-( 3.15[ ) (as in [3]) 



\\dMml + \\dMml<-{h\\\dMt)\\l 

1 



wdMmi) 



< 



m 



Ht)\ 



1/371, 



uit)\ 



<C{Eo + \\v{t)\\l + \\vit)\\l + \\u{t)\\ 

<C(£o + |Qo| + |Qo|'/' + Mo + Mi] 
1 



\dMt)\\ 



2h 



with C = C(a, /3,7) only depending on the parameters a,/3,7. Therefore 

\\dMt)\\l + WdMml < ^C{£o + |Qo| + |Qo|'/' + Mo + M3 + M^) K„. 



Note that Ko is independent of AI. Finally, since from (3.15) 

< 

\\uo\\l + \\dMt)\ 



vmi<2\\v{t)UdMm2<\\vmi + \\dMtm 



1 



\a\ 



we obtain 
(3.16) 



Ht)\[ 



< 



1 



|Qo 



Mo+ (l + |7/«|)ifo)'^' 



■■=K, 



with K independent of AI but depending on the initial data. Therefore, if the 
truncation parameter AI in (3.2) satisfies M > i^, we conclude by (3.16) and the 
definition of /*^, that (u, v) := {u^^ , v^) is actually a solution of the Schrodin- 



ger-KdV system ( 1.1 ). This completes the proof of Theorem 2.1 
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4. Proof of Proposition 13.21 

First of all, we need the following lemma concerning the global existence of the 
solution of the discrete problem (3.3). Due to the fact that the problem (3.3 1 is 
truncated, we are also able to obtain the essential estimate (4.1 ). For simplicity, 
we will omit the superscript M. 

Lemma 4.1. Let a,(3,j eR be such that > 0. Fix J e N, L > 0, h ^ L/{J+1), 
and {uq,Vq) € Xj^c x Xj^r. Then, for each T > 0, there exists a unique solution 

{u\v'') e C([0,T];Xjx) X C([0,T];Xj,k) 

of the problem ( |3.3[ ). Moreover, the following estimate is valid, 

\\D+u'^it)\\l + \\D+v'^it)\\l + \\u'\t)\\l + \\v'^it)\\l 

<C{\\u%,\\v^h,\\D+u'ih,\\D+v'oh,T,M), 

for all t e [0, T] . 



(4.1) 



Proof Let Sh{t) 



Ghit) 



be the unitary groups generated by the 



discrete operators iA'* and —D^ in the Xj space. The problem (3.3) can be written 
in the usual semigroup framework, as an integral equation in the x Xj^r space: 



(4.2) 



where 



u\t) 
v\t) 



Sh{t)u'^o+ / Sh[t-s)Js{u'\s),v'^{s))ds^:^i{u\v''), 



Gh{t)^^+ / Gh(t~s)JK{u''(s),v''{s))ds^:<^>2{u\v''), 



JK{u\v^) = -Dof{v^)+^D,{g'{v^)\u'^\''). 
For i? > ||uo II2 + ll'i'o II2 a-iid T > we consider the product space -B/jc x BJ^^, with 

- G C([0,r];Xj,c) : lkllL~([o,T];X,) < R} 

and similarly for _B|J ^. 

By (1.7 1, it is now a simple matter to prove that there exists T > such that 
the map 

{u\v'^) e X Bl^ ^ <i>{u\v'^) {<i>,{u\v'^),<i>2{u\v'^)) 

is a strict contraction on the complete metric space -B^c ^ ^rr- Thus, by the 
Banach fixed-point theorem, we obtain a unique local in time solution {u^,v^) of 
the problem ([33| in the space C([0, T]; Xj,c) x G[[G,T]; Xj^). 

To obtain a global solution, we must estimate the ^^-norm of u^[t) and v'^t), 
for each fixed h. From the equation (3.3a), the conservation of the P-norm of 
u'^{t) follows easily by taking the scalar product with and summation by parts. 
Applying the same procedure to equation (3.3b I, we find 

,7 



(4.3) 



1 



J 
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a„„ , o /37,i,./i/,M,4 



From the definition of /, the conservation of the P-norm of w'' and the fact that, 
for h fixed, || • ||oo < C!{h)\\ ■ II2, we derive that 

\\v\ml<\\vo\\l+Cih) f\\vHs)\\lds. 

Jo 

The conclusion now follows from a Gronwall argument. 

It remains to prove the inequality (4.1). In addition to the conservation of the 
P-norm of we have the conservation of the discrete energy: 

(4.4) 

To prove this identity, we proceed in the same way as in the continuous case |12j : 
Take the scalar product in Xj of the equation (3.3a) with dtu'\ take the real part, 
and use the equation (3.3b) and the skew-adjoint properties of the operators Dq 
and D^. 

Now we return to (4.3) and observe that from f{vj) = f{vj) — /(O) — f'{dj)vj, 
\g'\<l and (Ol we find 



J 



J 



dt\\vY2<CiM){\\D+v 
< C{M){\\D+v 



h\\2 



h\\2 
2 



\D+u' 



h\\2 



.'l||2^ 



2 ^ ll"^ 112;- 

Integrating on (0, t) and using Gronwall's lemma, we obtain 

(4.5) \\v''{t)h<ai{t,M) + a2{t,M) f \\D+v\s)\\l + \\D+u'\s)\\lds 

Jo 

for some continuous functions ai,a2. On the other hand, from the conservation of 
the energy (4.4) and using (1.5) we get 

WD+u'^ism + \\D+v'^{s)\\l < C{4,v^) + C,\\D+u'^is)h 

^^■^^ + C2 ^ hg{v,)\u,f + C3 ^ hF{v,). 

3=1 j=i 

But now, the definition of f^^ allows us (roughly) to bound F*^(u) by C{M)v'^. 
This is essential in view of the desired estimate ( |4.1[ ), since these terms would 
otherwise be cubic. We have 

Y,h\F{v,)\^Y.^\F{v,)-Fm = Y.h\j{e,)v,\ 

<Y.^f{e,f + \\v'^\\l 

for some 6j between and Vj. Now, 

i=l \vj\<M^ + l \vj\>M^ + l 
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Recall the definition of the truneated functions in (3.1 ). For the first sum, we have 
/(^j)^ < Vj < (M^ + ^y^'ji and for the second sum we have f{OjY < Thus we 
obtain 

J 



\2- 



Y^h\F{v,)\<C{M)\\v' 
Similarly, since the ^^-norm of m'' is conserved, we find 

Y,h\g{vMu,?<\\9\U\u'^\\l<C{M)\\u'^,\ 



h\\2 
2- 



J = l 



These estimates together with (4.5) and (4.6) give us 

\\D+^,\t)\\l + \\D+v\t)\\l < C{ulvlM) + C{M)\\v\t)\\l 

<C(ulvlM) + C{t,M) I \\D+uHs)\\l + \\D+v'\s)\\lds, 



where C{u^,v^,M) = C(||4'||2, P+M;f||2, Iko lb, ll^+^^o II2, M)- A Gronwall argu- 
ment, (4.5), and the conservation of ||u'*||2 give the conclusion (4.1 ). This completes 
the proof of Lemma |4.1| □ 



Proof of Proposition \3 .S\ We will use the interpolators P2,Pq defined in (1.2 1. 
Since d^Pi ~ P0-D+, it follows from the hypothesis P^Uq — >■ uq, P^Vq -> vq 
in H^{0,L), ([L4| and that 

\\P'lu\t)\\Hiio,L) < C, ||P^'^(t)|Ui(o,L) < C, te [0,T] 

with C = C(||uo||ffi, ll^'ollffi,?^,-^)- Thus, using the compactness of the embedding 
of H^{Q,L) into L^{0,L), we obtain, as /i — ?> (for a subsequence). 



(4.7) 



in L-([0,r];i7i(0,L)), 
P^u^^u, P'lv'^^v in L°"{[0,T];L\0,L)), 
for some u,v £ Hq{0,L). Also, we have from ( |1.3| ), 

(4.8) P'^u^^u, Pjfw'* -> w in L°°([0,r];L2(o,i)). 

To prove that u, v are solutions to the system ( |3.2p , we apply the piecewise constant 
interpolator Pq to the equations (3.3a|,(3.3b|: 

(4.9) idtP'^u'' + P'^A^u'' = /3P(j(|w''pu'') + aP'${g'^' {v^)u''') 

(4.10) dtP'^v'^ + P'^D'v'^ + P'^Dof^'iv'^) = 7PSi?o((5*')'(«')l"'n- 



From (1.3), (4.1) we have 

(4.11) Pifiv'') - P(,7(«") ^ in L-([0,T];L2(o,L)) 

and, since the piecewise constant interpolator commutes with nonliearities, it fol- 
lows from (4.8 1 that 

(4.12) P'ifiv'^) = fiPy ) ^ f{v) in L-([0,r];i2(o,L)). 
On the other hand, using ( |4.1[ ), 

Wd^P^lfiv'^m = |lPo"i^+/(«'')||2 = \\D+f iv'^)h < C{M), 
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and, from (4.11 ), (4.12 ) we deduce 



and so, 



P?/(«') ^ fiv) in L^{[0,T];H\0,L)) 
a,Pj/(«'') - in L-([0,r];i2(o,L)). 



Similarly, and using also (4.8), we find 

Pjfdu'^pw'') = |P(ju''|2pJfw'' ^ \ufu in L°°([0,r];L2(o,L)), 

P^giv'^'') = .9(P^'')PS"'' ^ 9{v)u in {[0,T]; L'{0, L)) , 

P^,D„{g\v'^)\u'^\') = d^P'l{g\v'^)\u'^\') - d^{g' {v)\u\^) in {[0,T]; L^{0, L)) , 

which allows us to pass to the limit on the corresponding terms in the weak formu- 
lation of the equations ( 1.1 ). 

It remains to analyze the terms PqA'^u'' and PqD^v'^. Let (j) £ 2^(0, be a 
test function. We have 



.7-1 



4>{x) dx 



= 22 J ^ {^(^ - - -h) + (j){x) - (pix + h)) dx 



3=2 

and so, by Taylor expansion of 0, 

,7-1 



J-1 



1/2 



i=2 



i=2 



Hence, from (4.1) we obtain 

(4.13) \\P^D\'^\\L^^[o,nH-Ho.L)) < C. 

If we now take a test function (p E D{(0,T) x (0, L)), we may compute in the sense 
of distributions 



{Pl^D^v'',ip) = ^D^Vj ip{x,t)dxdt 

Jo Jxj 



T J 

/ {dtif + 0{h)) dx dt 



as /i — > 0. Hence, we deduce from (4.131 that 

P^D^yh A in ([0, T] ; H'^ (0, L)) . 
In a similar way we prove that 



ph^h^h A Am in L°° {[Q,T\; R-"^ {Q, L)) 
and using the equations. 



idtP^u'' ^ idtu in L°°{[Q,T];H~'^{Q,L)) 
dtP'^v'^ - a^z; in ([0, T] ; i7-3(o, L)) . 
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Therefore, taking the Umit ft, — >■ in the weak formulation of the equations (4.9 ),(4.10 ) 
we obtain a weak solution 

(L-([0,T];i/i(0,L)))' 

n (C([0,T];iJ-2(o,L)) X C{[0,T];H-\0,L))), T > 0, 

of the proble m (|3.2[ ). To prove (3.11 1, recall that this solution satisfies the integral 
system (3.6|,(3.7|. Since 

\\miu,v)\\H^ <Ci\\uo\\„i,\\vo\\„i), 

\Hu,v)\\hi < Ci\\uo\\mA\vo\\H^), 
we deduce from ( 4.13| ) that 

{u,v) £ C{[0,T];H') X C{[0,T];L^). 



It remains to prove the energy inequality (3.121. Let us write the discrete energy 
(4.4) in the form 



(4.14) 



-I 
2 ' 



^-r-,d^P'yml + ^\\P'^u'\t)\\i 



+ aj 



/ g{P'Sv'')\P!Su''\'^ dx - a F{P!Sv'')dx^ E''{0). 
Jo Jo 



Now we recall the weak and strong convergences (|4.7[),(4.8 1. From the last term on 



the left-hand side of ( |4.14| ), and since |F'(^)| < CjCr, we find 

\FiP'yit)) - Fivit))\dx < cwp'yit) - «(t)ii2ii(Po^"(t))' + v\t)\ 



and so 



F{P^^v''{t))dx 



F(v[t))dx. 



Note that here it is essential that the spatial domain (0, V) be bounded. Indeed, a 
version of the energy inequality p. 12 1 on the whole line cannot be obtained using 



the available convergences (4.7),( |4.8[ ), which are local in space. 

Finally, from the strong convergence in iJ^(0, L) of the initial data (uq , Uq) and 
using the lower semi-continuity of the norm, we easily obtain from (4.141, in 
the limit /i ^ 0, the conclusion ( |3.12| : £(i) < £(0), t G [0,r]. This completes the 
proof of Proposition |3.2[ □ 



5. Numerical experiments 

In this section, we present some numerical computations using a fully discrete 
version of the method (2.1 ). We emphasize that these simulations are for the sake 
of illustration of our convergence result only. In particular, it would be interesting 
to perform more extensive numerical tests, such as determining the order of con- 
vergence, or employing more sophisticated time discretizations, which we do not 
perform here. 
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Given some time step r > 0, a spatial mesh size h, and initial data (uoj , voj )j=o J+i i 



we consider for n > the following algorithm. Set u"^^^^ ^(u"^^ and solve 



(5.1) ^ 

- z;^) + + i^o(^;;.'+l)2 ^ i?o(Kn, 

for j = 0, . . . , J+1. This corresponds to a semi-implicit Crank-Nickolson scheme for 
the Schrodinger equation and a fully implicit Euler scheme for the KdV equation. 
Because of the nonlinear terms, we perform a Newton iteration at each time step 
with a tolerance of 10~^. At each Newton iteration we solve independently a 



tridiagonal system for the first equation of (5.1) by a standard direct method, and 
the pentadiagonal system issuing from the second equation by an LU decomposition 
method. 

5.1. Comparison with exact solutions. We now test our scheme and illustrate 
our convergence result. We will simulate the following system. 



(5.2) 



dtv + dlv + vdxv = ^dx\u[' 



which is the same as (1.1) for a special choice of the parameters, except for the 
quasilinear term in the KdV equation which is formally equivalent to ^dxiv^)- In 



[6] we can find the following exact traveling wave solutions to ( |5.2[ ), 
(u,w) = (e"^*e"'=/2^(a; - ct), V(a; - ct)), 



with 



V2c*(l + 6a) I2d 



cosh(vc*2/) cosh {vc*y) 



Here, a G [-1/6,0] and w e K are given, and 2c = 1 + ^1 + |(1 + 6a), c* = 
+ o;^. We chose a ~ —1/12 and ui = for the simulations below. This gives a 
traveling wave speed c — 0.996516. 

In Figure [l] we present (in logarithmic scale) the relative error computed at 
T = 5 for a time step t = 0.0001 as a function of the mesh size. The computational 
domain is [—20, 50] and the number of spatial points ranges from 500 to 2500. In 
Figure [2] we present the relative error for a similar computation with T — 30 
and a time step r = 0.0005. The number of spatial points ranges from 200 to 600. 

In can be seen that the exact solution is approximated rather well by our simple 
numerical scheme, especially bearing in mind that the simulations (which, again, 
serve only to illustrate our results) were performed in a few minutes on a laptop 
running at 2.4 GHz. 
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Figure 1. Relative l? error, T = 5, t = 0.0001, as a function of 
mesh size. 500 to 2500 spatial points. 
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